#LyX 2.0 created this file. For more info see http://www.lyx.org/
\lyxformat 413
\begin_document
\begin_header
\textclass article
\use_default_options true
\maintain_unincluded_children false
\language catalan
\language_package default
\inputencoding auto
\fontencoding global
\font_roman default
\font_sans default
\font_typewriter default
\font_default_family default
\use_non_tex_fonts false
\font_sc false
\font_osf false
\font_sf_scale 100
\font_tt_scale 100

\graphics default
\default_output_format default
\output_sync 0
\bibtex_command default
\index_command default
\float_placement H
\paperfontsize default
\spacing single
\use_hyperref true
\pdf_bookmarks true
\pdf_bookmarksnumbered false
\pdf_bookmarksopen false
\pdf_bookmarksopenlevel 1
\pdf_breaklinks false
\pdf_pdfborder false
\pdf_colorlinks false
\pdf_backref false
\pdf_pdfusetitle true
\papersize default
\use_geometry true
\use_amsmath 1
\use_esint 1
\use_mhchem 1
\use_mathdots 1
\cite_engine basic
\use_bibtopic false
\use_indices false
\paperorientation portrait
\suppress_date false
\use_refstyle 0
\index Índex
\shortcut idx
\color #008000
\end_index
\topmargin 2cm
\bottommargin 2cm
\secnumdepth 3
\tocdepth 3
\paragraph_separation indent
\paragraph_indentation default
\quotes_language english
\papercolumns 1
\papersides 1
\paperpagestyle default
\tracking_changes false
\output_changes false
\html_math_output 0
\html_css_as_file 0
\html_be_strict false
\end_header

\begin_body

\begin_layout Title
Guia de referència per l'anàlisi d'espectres amb IRAF
\begin_inset Newline newline
\end_inset

Universitat de Barcelona
\end_layout

\begin_layout Author
Subirà Vallès, Montserrat
\begin_inset Newline newline
\end_inset

Blanco Cuaresma, Sergi
\begin_inset Newline newline
\end_inset

Peña Pérez, Dídac
\begin_inset Newline newline
\end_inset

Ruiz Dern, Laura
\end_layout

\begin_layout Standard
\begin_inset Newpage pagebreak
\end_inset


\end_layout

\begin_layout Standard
\begin_inset CommandInset toc
LatexCommand tableofcontents

\end_inset


\end_layout

\begin_layout Standard
\begin_inset Newpage pagebreak
\end_inset


\end_layout

\begin_layout Section
Objectius
\end_layout

\begin_layout Standard
El present document és una recopilació de les accions portades a termes
 amb el programa IRAF per a la reducció d'imatges preses per a la pràctica
 d'espectres, corresponents a l'assignatura d'Astronomia observacional de
 la Universitat de Barcelona.
\end_layout

\begin_layout Standard
No obstant, el document ha sigut preparat per tal de que també pugui ser
 utilitzat com una guia de referència general per a la reducció d'imatges
 amb IRAF.
 Els espectres analitzats són d'objectes puntuals (estrelles).
 En cas de disposar d'espectres extensos com nebuloses, es recomana complementar
 la informació presentada amb altres documents de referència com 
\emph on

\begin_inset Quotes eld
\end_inset

A user's guide to reducing slit spectra with IRAF
\begin_inset Quotes erd
\end_inset


\emph default
 de Phil Massey et al (1992).
\end_layout

\begin_layout Section
\begin_inset CommandInset label
LatexCommand label
name "sec:Instrumentació"

\end_inset

Instrumentació
\end_layout

\begin_layout Standard
Les imatges, reduïdes i analitzades amb IRAF per a la redacció d'aquest
 document, han estat obtingudes mitjançant una CCD ST-8 amb les següents
 característiques:
\end_layout

\begin_layout Itemize
Xip: KODAK KAF-1603ME de 1530x1020 píxels, cadascun d’ells de 9x9
\begin_inset Formula $\mu m$
\end_inset


\end_layout

\begin_layout Itemize
Sensibilitat (guany): 2.5 electrons per ADU (
\emph on
Analogic Digital Unit
\emph default
, és a dir, cada unitat en la imatge digital ha estat produïda per 2.5 electrons)
 
\end_layout

\begin_layout Itemize

\emph on
Read-out noise
\emph default
 (RON): 15 electrons per píxel que equivalen a 
\begin_inset Formula $15electrons\frac{1ADU}{2.5electrons}=6ADU$
\end_inset

 per píxel.
\end_layout

\begin_layout Itemize
Bits per píxel: 16 bits
\end_layout

\begin_layout Itemize
Nombre màxim de comptes: 
\begin_inset Formula $2^{16}=65516$
\end_inset

 
\end_layout

\begin_layout Itemize
Temperatures de treball: 20-25 graus per sota de la temperatura ambient
 (màxim -25) 
\end_layout

\begin_layout Section
Preparació de l'entorn IRAF
\end_layout

\begin_layout Standard
IRAF requereix un entorn Unix per a la seva instal·lació.
 En aquest apartat es descriu la instal·lació en un sistema Ubuntu GNU/Linux.
\end_layout

\begin_layout Subsection
Instal·lació
\end_layout

\begin_layout Standard
Com usuari administrador en un sistema Ubuntu GNU/Linux, es requereix disposar
 de la shell 'csh' i crear l'usuari 'iraf':
\end_layout

\begin_layout Standard
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

apt-get install csh
\end_layout

\begin_layout Plain Layout

mkdir -p /iraf/iraf
\end_layout

\begin_layout Plain Layout

useradd -d /iraf/iraf/local -m -s /bin/csh iraf
\end_layout

\end_inset


\end_layout

\begin_layout Standard
Un cop creat l'usuari, descarregar el fitxer 'iraf.lnux.x86_64.gz'
\begin_inset Foot
status open

\begin_layout Plain Layout
http://iraf.noao.edu/
\end_layout

\end_inset

 (o 'iraf.lnux.x86.gz' si el sistema és de 32 bits) i guardar-ho a '/iraf/'
 per tal de fer la descompressió:
\end_layout

\begin_layout Standard
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

cd /iraf/iraf/
\end_layout

\begin_layout Plain Layout

cat /iraf/iraf.lnux.x86_64.gz | zcat | tar -xpf - 
\end_layout

\begin_layout Plain Layout

chown -R iraf:iraf /iraf 
\end_layout

\end_inset


\end_layout

\begin_layout Standard
A continuació s'edita el fitxer '/iraf/iraf/unix/hlib/install' per tal d'afegir,
 al principi de tot:
\end_layout

\begin_layout Standard
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

set nonomatch
\end_layout

\end_inset


\end_layout

\begin_layout Standard
Seguidament, es prepara l'entorn per a la instal·lació:
\end_layout

\begin_layout Standard
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

cd /iraf/iraf/unix/hlib
\end_layout

\begin_layout Plain Layout

csh
\end_layout

\begin_layout Plain Layout

setenv iraf /iraf/iraf/
\end_layout

\begin_layout Plain Layout

source $iraf/unix/hlib/irafuser.csh
\end_layout

\end_inset


\end_layout

\begin_layout Standard
Es recomanable primer fer una simulació d'instal·lació per tal de detectar
 si existeix algun tipus de problema:
\end_layout

\begin_layout Standard
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

./install -n
\end_layout

\end_inset


\end_layout

\begin_layout Standard
Si no ha aparegut cap missatge d'error, es deixen tots els paràmetres per
 defecte, i es pot procedir a la instal·lació real:
\end_layout

\begin_layout Standard
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

./install
\end_layout

\end_inset


\end_layout

\begin_layout Standard
En aquesta ocasió, es deixaran tots els paràmetres per defecte a excepció
 de:
\end_layout

\begin_layout Standard
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

Configure IRAF Networking on this machine?  (yes): no
\end_layout

\begin_layout Plain Layout

...
\end_layout

\begin_layout Plain Layout

Create a default tapecap file?  (yes): no
\end_layout

\end_inset


\end_layout

\begin_layout Standard
Un cop finalitzada la instal·lació d'IRAF, es surt de la shell 'csh':
\end_layout

\begin_layout Standard
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

% exit
\end_layout

\end_inset


\end_layout

\begin_layout Standard
Per a la versió 2.15.1a 64 bits, es recomana descarregar unes correccions
 de bugs que seran utilitzades per a alguna tasca concreta de fotometria:
\end_layout

\begin_layout Standard
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

cd /iraf/iraf/bin.linux64/
\end_layout

\begin_layout Plain Layout

mv ecl.e ecl64.e
\end_layout

\begin_layout Plain Layout

wget ftp://iraf.noao.edu/iraf/v215/support/linux/ecl.e
\end_layout

\begin_layout Plain Layout

mv ecl.e ecl32.e
\end_layout

\begin_layout Plain Layout

ln -s ecl64.e ecl.e 
\end_layout

\begin_layout Plain Layout

\end_layout

\begin_layout Plain Layout

cd /iraf/iraf/noao/bin.linux64
\end_layout

\begin_layout Plain Layout

mv x_photcal.e x_photcal64.e
\end_layout

\begin_layout Plain Layout

wget ftp://iraf.noao.edu/iraf/v215/support/linux/x_photcal.e
\end_layout

\begin_layout Plain Layout

chmod 755 x_photcal.e
\end_layout

\begin_layout Plain Layout

mv x_photcal.e x_photcal32.e
\end_layout

\begin_layout Plain Layout

ln -s x_photcal64.e x_photcal.e 
\end_layout

\begin_layout Plain Layout

\end_layout

\begin_layout Plain Layout

chown iraf:iraf /iraf/iraf/bin.linux64/* /iraf/iraf/noao/bin.linux64/*
\end_layout

\end_inset


\end_layout

\begin_layout Standard
A continuació cal instal·lar 'xgterm', la terminal que utilitza IRAF:
\end_layout

\begin_layout Standard
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

mkdir x11iraf-v2.0BETA-bin.linux/
\end_layout

\begin_layout Plain Layout

cd x11iraf-v2.0BETA-bin.linux/
\end_layout

\begin_layout Plain Layout

wget ftp://iraf.noao.edu/iraf/x11iraf/x11iraf-v2.0BETA-bin.linux.tar.gz
\end_layout

\begin_layout Plain Layout

tar -zxvf x11iraf-v2.0BETA-bin.linux.tar.gz
\end_layout

\begin_layout Plain Layout

sudo ./install
\end_layout

\end_inset


\end_layout

\begin_layout Standard
Un cop instal·lat, serà necessari crear el fitxer '/usr/local/bin/iraf.sh'
 amb el següent contingut:
\end_layout

\begin_layout Standard
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

#!/bin/bash
\end_layout

\begin_layout Plain Layout

cd $HOME
\end_layout

\begin_layout Plain Layout

/usr/local/bin/ds9 &
\end_layout

\begin_layout Plain Layout

/usr/local/bin/xgterm -G 1100x350+10-10 -sbr -fn 9x15 -sl 400 -cr red -title
 IRAF -bg "black" -fg "white" -geometry 120x34 -e cl
\end_layout

\end_inset


\end_layout

\begin_layout Standard
I donar-li permisos d'execució:
\end_layout

\begin_layout Standard
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

chmod 755 /usr/local/bin/iraf.sh
\end_layout

\end_inset


\end_layout

\begin_layout Standard
Finalment, es descarrega el programa DS9 (entorn per a la visualització
 i anàlisi d'imatges FITS) 
\begin_inset Foot
status open

\begin_layout Plain Layout
http://hea-www.harvard.edu/RD/ds9/
\end_layout

\end_inset

 per tal de descomprimir-lo i copiar l'executable a un dels directoris del
 sistema:
\end_layout

\begin_layout Standard
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

wget http://hea-www.harvard.edu/saord/download/ds9/linux64/ds9.linux64.6.2.tar.gz
\end_layout

\begin_layout Plain Layout

tar -zxvf ds9.linux64.6.2.tar.gz
\end_layout

\begin_layout Plain Layout

mv ds9 /usr/local/bin
\end_layout

\end_inset


\end_layout

\begin_layout Subsection
Execució
\end_layout

\begin_layout Standard
Per poder accedir a IRAF amb un usuari normal, la primera vegada s'haurà
 d'executar 'mkiraf' i seleccionar 'xgterm' com a terminal:
\end_layout

\begin_layout Standard
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

cd && mkiraf && iraf.sh
\end_layout

\end_inset


\end_layout

\begin_layout Standard
D'aquesta forma 'mkiraf' crearà el fitxer 'login.cl' i el directori 'uparm',
 necessaris per IRAF.
 Les següents vegades que es vulgui accedir a IRAF només serà necessari
 executar:
\end_layout

\begin_layout Standard
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

iraf.sh
\end_layout

\end_inset


\end_layout

\begin_layout Subsection
Comandes bàsiques
\end_layout

\begin_layout Standard
Algunes comandes d'utilitat general que poden ser utilitzades en iniciar
 l'entorn:
\end_layout

\begin_layout Description
logout Sortir
\end_layout

\begin_layout Description
package Llista paquets que poden ser carregats escrivint directament el
 seu nom (p.ex.
 imred, ccdred)
\end_layout

\begin_layout Description
bye Descarrega el paquet actual
\end_layout

\begin_layout Description
help Ajuda de comandes concretes
\end_layout

\begin_layout Description
show
\begin_inset space ~
\end_inset

|
\begin_inset space ~
\end_inset

page Llista les variables d'entorn
\end_layout

\begin_layout Description
show
\begin_inset space ~
\end_inset

editor Mostra la variable d'entorn 'editor'
\end_layout

\begin_layout Description
set
\begin_inset space ~
\end_inset

editor=vim Canvia una variable d'entorn (en aquest cas fa referència a l'editor
 que s'executarà quan s'utilitzi la comanda 'edit fitxer')
\end_layout

\begin_layout Description
cd
\begin_inset space ~
\end_inset

imdir Canvia al directori d'imatges FIT de l'usuari (habitualment /iraf/imdir/
\emph on
usuari
\emph default
/)
\end_layout

\begin_layout Description
pwd Mostra el directori actual
\end_layout

\begin_layout Description
ls Llista fitxers i directoris
\end_layout

\begin_layout Description
lparam
\begin_inset space ~
\end_inset

imstat Llista els paràmetres de la comanda 'imstat'.
\end_layout

\begin_layout Description
epar
\begin_inset space ~
\end_inset

imstat Modifica els paràmetres de la comanda 'imstat' que pot ser executada
 a continuació.
 Per sortir de l'edició guardant CTRL+D, sense guardar CTRL+C.
\end_layout

\begin_layout Description
imstat
\begin_inset space ~
\end_inset

fitxer.fits Mostra estadístiques del fitxer FIT.
 És important que aquest tingui extensió '.fits' o '.fits' en minúscules.
\end_layout

\begin_layout Description
delete
\begin_inset space ~
\end_inset

fitxer Esborra un fitxer
\end_layout

\begin_layout Description
cl
\begin_inset space ~
\end_inset

<
\begin_inset space ~
\end_inset

calibrate.cl Executa un script anomenat 'calibrate.cl' ubicat al directori
 actual
\end_layout

\begin_layout Description
less
\begin_inset space ~
\end_inset

fitxer Visualitza un fitxer de text ('q' per sortir)
\end_layout

\begin_layout Description
!
\begin_inset space ~
\end_inset

gedit
\begin_inset space ~
\end_inset

fitxer Executa la comanda del sistema 'less fitxer'
\end_layout

\begin_layout Section
Reducció d'imatges
\end_layout

\begin_layout Standard
A continuació es descriu el procés de reducció d'imatges brutes per tal
 de poder extreure els espectres.
\end_layout

\begin_layout Standard
IRAF requereix treballar amb imatges FITS, en cas de tenir captures en format
 ST8 serà necessari fer una conversió, per exemple, fent ús del programa
 FITSConverter
\begin_inset Foot
status open

\begin_layout Plain Layout
http://www.sbig.com/sbwhtmls/softpage.htm
\end_layout

\end_inset

:
\end_layout

\begin_layout Standard
\begin_inset Float figure
wide false
sideways false
status open

\begin_layout Plain Layout
\align center
\begin_inset Graphics
	filename Images/FITSConverter.png
	lyxscale 50
	scale 35

\end_inset


\end_layout

\begin_layout Plain Layout
\begin_inset Caption

\begin_layout Plain Layout
FITSConverter, només disponible per sistemes Microsoft Windows.
\end_layout

\end_inset


\end_layout

\end_inset


\end_layout

\begin_layout Standard
Cal destacar que les imatges FITS han de tenir extensions 'fit' o 'fits'
 en minúscules ja que sinó és probable que l'IRAF no les pugui obrir correctamen
t.
 Per renombrar massivament un conjunt de fitxers a GNU/Linux s'empra:
\end_layout

\begin_layout Standard
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

$ rename 's/
\backslash
.FIT$/.fits/' *
\end_layout

\begin_layout Plain Layout

$ rename 's/
\backslash
.FITS$/.fits/' *
\end_layout

\end_inset


\end_layout

\begin_layout Subsection
Preparació dels fitxers d'entrada
\end_layout

\begin_layout Enumerate
Abans de començar, és important guardar una còpia de totes les imatges FITS
 que s'utilitzaran donat que el procés sobreescriurà moltes d'elles.
 
\end_layout

\begin_layout Enumerate
S'inicia el programa gràfic DS9 mitjançant l'execució de 'ds9':
\begin_inset Float figure
wide false
sideways false
status open

\begin_layout Plain Layout
\align center
\begin_inset Graphics
	filename Images/ds9.png
	lyxscale 50
	scale 35

\end_inset


\end_layout

\begin_layout Plain Layout
\begin_inset Caption

\begin_layout Plain Layout
DS9, visualitzador d'imatges FITS
\end_layout

\end_inset


\end_layout

\begin_layout Plain Layout

\end_layout

\end_inset


\end_layout

\begin_layout Enumerate
S'inicia l'entorn IRAF, situant-se al directori a on estan les imatges (habitual
ment 
\begin_inset Quotes eld
\end_inset

/iraf/imdir/
\emph on
usuari
\emph default
/
\begin_inset Quotes erd
\end_inset

, referenciat per la variable 'imdir'):
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

$ cl
\end_layout

\begin_layout Plain Layout

ecl>
\end_layout

\begin_layout Plain Layout

ecl> cd imdir
\end_layout

\end_inset


\begin_inset Float figure
wide false
sideways false
status open

\begin_layout Plain Layout
\align center
\begin_inset Graphics
	filename Images/iraf.png
	lyxscale 50
	scale 35

\end_inset


\end_layout

\begin_layout Plain Layout
\begin_inset Caption

\begin_layout Plain Layout
Entorn IRAF
\end_layout

\end_inset


\end_layout

\begin_layout Plain Layout

\end_layout

\end_inset

Pel present treball es disposa dels següents fitxers inicials (entre altres):
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

ecl> ls
\end_layout

\begin_layout Plain Layout

database			      grup4_alta_obj_HIP28634_90_1.fits
\end_layout

\begin_layout Plain Layout

grup1_alta_lamp_1_2.fits	      grup4_alta_obj_HIP28634_90_2.fits
\end_layout

\begin_layout Plain Layout

grup1_alta_obj_HD81192_30_1.fits      grup4_alta_obj_HIP28634_90_3.fits
\end_layout

\begin_layout Plain Layout

grup1_alta_obj_hip57939_60_1.fits     grup4_alta_obj_HIP35136_60_1.fits
\end_layout

\begin_layout Plain Layout

grup1_alta_obj_SAO77800_60_1.fits     grup4_baixa_lamp_10_2.fits
\end_layout

\begin_layout Plain Layout

grup1_alta_obj_SAO77800_60_2.fits     grup4_baixa_obj_Arcturus_10_1.fits
\end_layout

\begin_layout Plain Layout

grup1_alta_obj_SAO77800_60_3.fits     grup4_baixa_obj_Denebola_10_1.fits
\end_layout

\begin_layout Plain Layout

grup1_baixa_lamp_1_2.fits	      grup4_baixa_obj_Nekkar_10_1.fits
\end_layout

\begin_layout Plain Layout

grup1_dark_30_1.fits		      grup4_dark_10_1.fits
\end_layout

\begin_layout Plain Layout

grup1_dark_30_2.fits		      grup4_dark_2_1.fits
\end_layout

\begin_layout Plain Layout

grup1_dark_30_combined.fits	      grup4_dark_60_1.fits
\end_layout

\begin_layout Plain Layout

grup2_alta_lamp_10_1.fits	      grup4_dark_90_1.fits
\end_layout

\begin_layout Plain Layout

grup2_alta_obj_SAO58636_30_1.fits     grup5_alta_lamp_10_2.fits
\end_layout

\begin_layout Plain Layout

grup2_alta_obj_SAO60896_150_1.fits    grup5_baixa_lamp_10_2.fits
\end_layout

\begin_layout Plain Layout

grup2_alta_obj_SAO60896_150_2.fits    grup5_baixa_obj_alphard_10_1.fits
\end_layout

\begin_layout Plain Layout

grup3_dark_30_combined.fits	      scripts
\end_layout

\begin_layout Plain Layout

...
\end_layout

\end_inset

Totes les imatges presenten un prefix 
\begin_inset Quotes eld
\end_inset

grup
\begin_inset Quotes erd
\end_inset

 per tal d'identificar quin equip de treball les ha obtingut.
 A continuació, s'especifica si es tracta d'una imatge d'alta o baixa resolució,
 seguint del text 'dark', 'obj' o 'lamp' en funció de si es un dark, una
 estrella o una imatge de lampara.
 Finalment, s'especifica el nomb de l'estrella pels casos que aplica i els
 segons d'exposició.
\end_layout

\begin_layout Enumerate
Per visualitzar les imatges al programa DS9, es pot utilitzar la comanda
 'display':
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

ecl> display grup3_baixa_obj_betelgeuse_30_1.fits 1 fill+
\end_layout

\end_inset

DS9 disposa de diferents 'frames' a on es poden carregar diferents imatges
 simultàniament i alternar-les entre elles pressionant 'tabulador':
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

ecl> display grup3_baixa_obj_betelgeuse_30_1.fits 1 fill+
\end_layout

\begin_layout Plain Layout

ecl> display grup3_baixa_obj_castor_30_1.fits 2 fill+
\end_layout

\begin_layout Plain Layout

ecl> display grup3_baixa_obj_meissa_30_2.fits 3 fill+
\end_layout

\end_inset

En cas de que la visualització no sigui satisfactòria, es poden ajustar
 diferents paràmetres al menú 'Scale' de DS9.
\end_layout

\begin_layout Enumerate
Per defecte, les funcions d'extracció d'espectres esperen que l'eix de dispersió
 en posició vertical.
 En aquest sentit es poden rotar totes les imatges:
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

ecl> imtranspose *.fits *.fits
\end_layout

\end_inset


\begin_inset Float figure
wide false
sideways false
status open

\begin_layout Plain Layout
\align center
\begin_inset Graphics
	filename Images/Spectra-horitzontal.png
	lyxscale 70
	scale 33

\end_inset

 
\begin_inset Graphics
	filename Images/Spectra-Vertical.png
	lyxscale 50
	scale 25

\end_inset


\end_layout

\begin_layout Plain Layout
\begin_inset Caption

\begin_layout Plain Layout
Transposició d'una imatge per a que l'espectre (eix de dispersió) sigui
 vertical.
\end_layout

\end_inset


\end_layout

\end_inset

Una alternativa més eficient és indicar a la capçalera dels fitxers l'orientació
 de l'eix de dispersió (1 - horitzontal, 2 - vertical):
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

hedit *.fits DISPAXIS 1 add+
\end_layout

\end_inset


\end_layout

\begin_layout Enumerate
Per tal de llistar la capçalera d'una de les imatges a reduir s'utilitzarà
 'imhead':
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

ecl> imhead grup3_baixa_obj_betelgeuse_30_1.fits l+
\end_layout

\begin_layout Plain Layout

\end_layout

\begin_layout Plain Layout

grup3_baixa_obj_betelgeuse_30_1.fits[1530,1020][ushort]: 
\end_layout

\begin_layout Plain Layout

No bad pixels, min=0., max=65535.
\end_layout

\begin_layout Plain Layout

Line storage mode, physdim [1530,1020], length of user area 1580 s.u.
\end_layout

\begin_layout Plain Layout

Created Thu 02:00:00 02-Jun-2011, Last modified Thu 20:52:07 02-Jun-2011
\end_layout

\begin_layout Plain Layout

Pixel file "grup3_baixa_obj_betelgeuse_30_1.fits" [ok]
\end_layout

\begin_layout Plain Layout

EXTEND  =                    T / FITS dataset may contain extensions
\end_layout

\begin_layout Plain Layout

COMMENT = 'SBIG FITS header format per:'
\end_layout

\begin_layout Plain Layout

COMMENT = ' http://www.sbig.com/pdffiles/SBFITSEXT_1r0.pdf'
\end_layout

\begin_layout Plain Layout

BZERO   =                32768 / offset data range to that of unsigned short
\end_layout

\begin_layout Plain Layout

BSCALE  =                    1 / default scaling factor
\end_layout

\begin_layout Plain Layout

OBJECT  = '        '
\end_layout

\begin_layout Plain Layout

TELESCOP= '        '
\end_layout

\begin_layout Plain Layout

INSTRUME= 'SBIG ST-8 Dual CCD Camera' / Camera Model
\end_layout

\begin_layout Plain Layout

OBSERVER= 'Xavier Franch'
\end_layout

\begin_layout Plain Layout

DATE-OBS= '2011-04-05T21.00.32.000' / GMT START OF EXPOSURE
\end_layout

\begin_layout Plain Layout

EXPTIME =                  30.
 / EXPOSURE IN SECONDS
\end_layout

\begin_layout Plain Layout

CCD-TEMP=                 0.19 / CCD TEMP IN DEGREES C
\end_layout

\begin_layout Plain Layout

XPIXSZ  =                   9.
 / PIXEL WIDTH IN MICRONS
\end_layout

\begin_layout Plain Layout

YPIXSZ  =                   9.
 / PIXEL HEIGHT IN MICRONS
\end_layout

\begin_layout Plain Layout

XBINNING=                    1 / HORIZONTAL BINNING FACTOR
\end_layout

\begin_layout Plain Layout

YBINNING=                    1 / VERTICAL BINNING FACTOR
\end_layout

\begin_layout Plain Layout

XORGSUBF=                    0 / SUB_FRAME ORIGIN X_POS
\end_layout

\begin_layout Plain Layout

YORGSUBF=                    0 / SUB_FRAME ORIGIN Y_POS
\end_layout

\begin_layout Plain Layout

EGAIN   =                 2.57 / ELECTRONS PER ADU
\end_layout

\begin_layout Plain Layout

FOCALLEN=            4045.9914 / FOCAL LENGTH IN MM
\end_layout

\begin_layout Plain Layout

APTDIA  =                101.6 / APERTURE DIAMETER IN MM
\end_layout

\begin_layout Plain Layout

APTAREA =        125200.007664 / APERTURE AREA IN SQ-MM
\end_layout

\begin_layout Plain Layout

CBLACK  =                  115 / BLACK ADU FOR DISPLAY
\end_layout

\begin_layout Plain Layout

CWHITE  =                  365 / WHITE ADU FOR DISPLAY
\end_layout

\begin_layout Plain Layout

PEDESTAL=                 -100 / ADD TO ADU FOR 0-BASE
\end_layout

\begin_layout Plain Layout

DATAMAX =                65535 / SATURATION LEVEL
\end_layout

\begin_layout Plain Layout

SBSTDVER= 'SBFITSEXT Version 1.0' / SBIG FITS EXTENSIONS VER
\end_layout

\begin_layout Plain Layout

SWACQUIR= 'WinOPS Ver 5.41 Build 6-NT' / DATA ACQ SOFTWARE
\end_layout

\begin_layout Plain Layout

SWCREATE= 'WinOPS Ver 5.41 Build 6-NT'
\end_layout

\begin_layout Plain Layout

FILTER  = 'Red     '           / OPTICAL FILTER NAME
\end_layout

\begin_layout Plain Layout

SNAPSHOT=                    1 / NUMBER IMAGES COADDED
\end_layout

\begin_layout Plain Layout

DATE    = '2011-06-02'         / GMT DATE WHEN THIS FILE CREATED
\end_layout

\begin_layout Plain Layout

RESMODE =                    6 / RESOLUTION MODE
\end_layout

\begin_layout Plain Layout

EXPSTATE= '125     '           / EXPOSURE STATE (HEX)
\end_layout

\begin_layout Plain Layout

RESPONSE=                3000.
 / CCD RESPONSE FACTOR
\end_layout

\begin_layout Plain Layout

NOTE    = 'Local time:4/5/2011 at 21:00:32'
\end_layout

\begin_layout Plain Layout

TRAKTIME=                   0.
 / TRACKING EXPOSURE
\end_layout

\begin_layout Plain Layout

SWMODIFY= 'SBIG Batch SBIG to FITS Converter'
\end_layout

\begin_layout Plain Layout

HISTORY Convert from SBIG to FITS Format
\end_layout

\begin_layout Plain Layout

\end_layout

\end_inset

S'observa que la imatge és de 1530x1020 píxels i que no existeixen les següents
 seccions:
\end_layout

\begin_deeper
\begin_layout Enumerate
TRIMSEC: Píxels útils de la imatge
\end_layout

\begin_layout Enumerate
BIASEC: Regió no exposada a la llum (
\emph on
overscan
\emph default
)
\end_layout

\end_deeper
\begin_layout Enumerate
Es procedeix a analitzar la imatge localitzant la coordenada 'y' a on es
 troba la part central de l'espectre, per després executar 'implot' i selecciona
r la línia 'y' (p.ex.
 teclejant ':l 390'):
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

display grup3_baixa_obj_betelgeuse_30_1.fits 1 fill+
\end_layout

\begin_layout Plain Layout

implot grup3_baixa_obj_betelgeuse_30_1.fits
\end_layout

\end_inset


\begin_inset Float figure
wide false
sideways false
status open

\begin_layout Plain Layout
\align center
\begin_inset Graphics
	filename Images/implot-betelgueuse-390-without-dark.png
	lyxscale 50
	scale 25

\end_inset


\begin_inset Graphics
	filename Images/implot-betelgueuse-390-with-dark.png
	lyxscale 50
	scale 25

\end_inset


\end_layout

\begin_layout Plain Layout
\begin_inset Caption

\begin_layout Plain Layout
Gràfics centrats a la línia 390 de les imatges de baixa resolució de Betelgeuse
 original (esquerra) i amb darks aplicats (dreta).
\end_layout

\end_inset


\end_layout

\begin_layout Plain Layout

\end_layout

\end_inset

 Les comandes que es poden utilitzar directament sobre el gràfic són:
\end_layout

\begin_deeper
\begin_layout Enumerate
\begin_inset Quotes eld
\end_inset

:l 300
\begin_inset Quotes erd
\end_inset

: visualitza els valors de les columnes per a la línia 300
\end_layout

\begin_layout Enumerate
\begin_inset Quotes eld
\end_inset

:c 300
\begin_inset Quotes erd
\end_inset

: visualitza els valors de les línies per a la columna 300
\end_layout

\begin_layout Enumerate
\begin_inset Quotes eld
\end_inset

C
\begin_inset Quotes erd
\end_inset

: posició del cursor
\end_layout

\begin_layout Enumerate
\begin_inset Quotes eld
\end_inset

ee
\begin_inset Quotes erd
\end_inset

: amplia la regió a on està situat el cursor
\end_layout

\begin_layout Enumerate
\begin_inset Quotes eld
\end_inset

aa
\begin_inset Quotes erd
\end_inset

: torna a mostrar la gràfica a escala
\end_layout

\begin_layout Enumerate
\begin_inset Quotes eld
\end_inset

q
\begin_inset Quotes erd
\end_inset

: sortir
\end_layout

\end_deeper
\begin_layout Enumerate
Es defineixen els camps de capçalera que s'han detectat que faltaven:
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

hedit *.fits TRIMSEC [1:1530,1:1020] verify- show+ add+
\end_layout

\begin_layout Plain Layout

hedit *.fits BIASEC "[]" verify- show+ add+
\end_layout

\end_inset


\end_layout

\begin_layout Enumerate
Els darks, flats i imatges originals han de disposar d'un camp IMAGETYP
 que les identifiqui:
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

hedit *_dark_*.fits IMAGETYP dark verify- show+ add+
\end_layout

\begin_layout Plain Layout

hedit grup*_alta_*.fits IMAGETYP object verify- show+ add+
\end_layout

\begin_layout Plain Layout

hedit grup*_baixa_*.fits IMAGETYP object verify- show+ add+
\end_layout

\begin_layout Plain Layout

hedit grup*_lamp_*.fits IMAGETYP comp verify- show+ add+ 
\end_layout

\end_inset


\end_layout

\begin_layout Subsection
Reducció de les imatges
\end_layout

\begin_layout Standard
El procés per a la reducció d'una imatge requereix principalment portar
 a terme les següents operacions:
\end_layout

\begin_layout Standard
\align center
\begin_inset Formula 
\[
Imatge\, reduida=\frac{\left(Imatge\, original_{filtre}-Dark_{temps\, imatge}\right)}{\left(Flat_{filtre}-Dark_{temps\, flat}\right)_{normalitzat}}
\]

\end_inset


\end_layout

\begin_layout Standard
No obstant, per a l'obtenció d'un espectre serà suficient amb restar els
 darks i no considerar l'eficiència individual dels píxels (flats):
\end_layout

\begin_layout Standard
\align center
\begin_inset Formula 
\[
Imatge\, reduida=Imatge\, original_{filtre}-Dark_{temps\, imatge}
\]

\end_inset


\end_layout

\begin_layout Enumerate
Per tal de que la tasca de reducció 'ccdproc' interpreti correctament les
 capçaleres dels fitxers, es configura per realitzar fotometria amb CCD
 directe:
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

ecl> imred
\end_layout

\begin_layout Plain Layout

imred> ccdred
\end_layout

\begin_layout Plain Layout

imred> cd imdir
\end_layout

\begin_layout Plain Layout

ccdred> setinstrument direct
\end_layout

\end_inset

Es modifica el paràmetre 'pixelty' de manera que les operacions es realitzin
 utilitzant coma flotant però les imatges resultants continuïn sent representade
s per nombres reals:
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

(pixelty=           short real) Output and calculation pixel datatypes 
\end_layout

\end_inset

La resta de paràmetres no es modifiquen (ctrl+d, ctrl+d).
\end_layout

\begin_layout Enumerate
Els fitxers dark de la mateixa exposició i presses a la mateixa nit es combinen
 mitjançant les següents comandes:
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

darkcombine grup1_dark_30_*.fits output=grup1_dark_30_combined.fits combine=averag
e reject=minmax ccdtype=dark scale=exposure process-
\end_layout

\begin_layout Plain Layout

darkcombine grup3_dark_10_*.fits output=grup3_dark_10_combined.fits combine=averag
e reject=minmax ccdtype=dark scale=exposure process-
\end_layout

\begin_layout Plain Layout

darkcombine grup3_dark_30_*.fits output=grup3_dark_30_combined.fits combine=averag
e reject=minmax ccdtype=dark scale=exposure process-
\end_layout

\begin_layout Plain Layout

darkcombine grup6_dark_1_*.fits output=grup6_dark_1_combined.fits combine=average
 reject=minmax ccdtype=dark scale=exposure process-
\end_layout

\begin_layout Plain Layout

darkcombine grup6_dark_10_*.fits output=grup6_dark_10_combined.fits combine=averag
e reject=minmax ccdtype=dark scale=exposure process-
\end_layout

\begin_layout Plain Layout

darkcombine grup6_dark_60_*.fits output=grup6_dark_60_combined.fits combine=averag
e reject=minmax ccdtype=dark scale=exposure process-
\end_layout

\end_inset

D'aquesta forma es redueixen possibles efectes de rajos còsmics que hagin
 caigut sobre la CCD al moment de l'exposició.
\end_layout

\begin_layout Enumerate
Finalment, es procedeix a realitzar la reducció de les imatges:
\end_layout

\begin_deeper
\begin_layout Standard
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

# Grup 2
\end_layout

\begin_layout Plain Layout

ccdproc grup2*_alta_*_10_*.fits ccdtype=object oversca- trim- zerocor- darkcor+
 flatcor- dark=grup2_dark_10_1.fits
\end_layout

\begin_layout Plain Layout

ccdproc grup2*_alta_*_30_*.fits ccdtype=object oversca- trim- zerocor- darkcor+
 flatcor- dark=grup2_dark_30_1.fits
\end_layout

\begin_layout Plain Layout

ccdproc grup2*_alta_*_60_*.fits ccdtype=object oversca- trim- zerocor- darkcor+
 flatcor- dark=grup2_dark_60_1.fits
\end_layout

\begin_layout Plain Layout

ccdproc grup2*_alta_*_120_*.fits ccdtype=object oversca- trim- zerocor- darkcor+
 flatcor- dark=grup2_dark_120_1.fits
\end_layout

\begin_layout Plain Layout

ccdproc grup2*_alta_*_150_*.fits ccdtype=object oversca- trim- zerocor- darkcor+
 flatcor- dark=grup2_dark_150_1.fits
\end_layout

\begin_layout Plain Layout

\end_layout

\begin_layout Plain Layout

ccdproc grup2*_baixa_*_10_*.fits ccdtype=object oversca- trim- zerocor- darkcor+
 flatcor- dark=grup2_dark_10_1.fits
\end_layout

\begin_layout Plain Layout

ccdproc grup2*_baixa_*_15_*.fits ccdtype=object oversca- trim- zerocor- darkcor+
 flatcor- dark=grup2_dark_15_1.fits
\end_layout

\begin_layout Plain Layout

ccdproc grup2*_baixa_*_30_*.fits ccdtype=object oversca- trim- zerocor- darkcor+
 flatcor- dark=grup2_dark_30_1.fits
\end_layout

\begin_layout Plain Layout

ccdproc grup2*_baixa_*_45_*.fits ccdtype=object oversca- trim- zerocor- darkcor+
 flatcor- dark=grup2_dark_45_1.fits
\end_layout

\begin_layout Plain Layout

\end_layout

\begin_layout Plain Layout

# Grup 3
\end_layout

\begin_layout Plain Layout

ccdproc grup3*_baixa_*_10_*.fits ccdtype=object oversca- trim- zerocor- darkcor+
 flatcor- dark=grup3_dark_10_combined.fits
\end_layout

\begin_layout Plain Layout

ccdproc grup3*_baixa_*_30_*.fits ccdtype=object oversca- trim- zerocor- darkcor+
 flatcor- dark=grup3_dark_30_combined.fits
\end_layout

\end_inset


\end_layout

\end_deeper
\begin_layout Section
Espectres
\end_layout

\begin_layout Standard
En aquesta secció s'exposen els passos necessaris per extreure i calibrar
 l'espectre d'una estrella a partir de les imatges reduïdes.
\end_layout

\begin_layout Subsection
\begin_inset CommandInset label
LatexCommand label
name "sub:Característiques-de-les"

\end_inset

Característiques de les imatges
\end_layout

\begin_layout Standard
A continuació s'indiquen algunes de les funcions disponibles per analitzar
 les imatges i obtenir valors d'interès per a l'extracció final de l'espectre.
\end_layout

\begin_layout Enumerate
Amb l'objectiu d'analitzar la traça de l'estrella, es visualitza la imatge
 corresponent utilitzant 'display' i 'implot':
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

display grup3_baixa_obj_betelgeuse_30_1.fits 1 fill+
\end_layout

\begin_layout Plain Layout

implot grup3_baixa_obj_betelgeuse_30_1.fits
\end_layout

\end_inset


\end_layout

\begin_layout Enumerate
S'indica que la gràfica d'implot mostri una línia que talli l'espectre (p.ex.
 ':l 500') i a continuació:
\end_layout

\begin_deeper
\begin_layout Enumerate
Es situa el cursor per sota de l'eix X abans del pic de flux (ubicació de
 l'espectre de l'estrella) i es pressiona 'p'.
\end_layout

\begin_layout Enumerate
Es desplaça el cursor per situar-ho després del pic i es torna a pressionar
 'p'.
 
\begin_inset Float figure
wide false
sideways false
status open

\begin_layout Plain Layout
\align center
\begin_inset Graphics
	filename Images/implot-betelgueuse-caracterize.png
	lyxscale 50
	scale 35

\end_inset


\end_layout

\begin_layout Plain Layout
\begin_inset Caption

\begin_layout Plain Layout
Gràfics centrats a la línia 390 de les imatges de baixa resolució de Betelgeuse
 original (esquerra) i amb darks aplicats (dreta).
\end_layout

\end_inset


\end_layout

\end_inset

D'aquesta forma s'obté el FWHM (
\emph on
Full width at half maximum
\emph default
), es a dir, l'amplada en píxels a alçada meitat.
 Si en comptes de 'p' s'utilitza 's', llavors es pot extreure el número
 de píxels entre els dos punts i altres valors estadístics.
\end_layout

\end_deeper
\begin_layout Enumerate
A continuació es valida el grau de perpendicularitat de l'espectre respecte
 a l'eix Y.
 S'escolleixen 3 linies a principi, meitat i final de l'espectre (coordenada
 Y visualitzada al ds9) i des de l'implot:
\end_layout

\begin_deeper
\begin_layout Enumerate
':l 450' Visualitza la línia a meitat de l'espectre (p.ex.
 450)
\end_layout

\begin_layout Enumerate
':o' No esborrar la línia anterior
\end_layout

\begin_layout Enumerate
':l 150' Visualitza la línia a l'inici de l'espectre (p.ex.
 150)
\end_layout

\begin_layout Enumerate
':o' No esborrar la línia anterior
\end_layout

\begin_layout Enumerate
':l 750' Visualitza la línia al final de l'espectre (p.ex.
 750)
\begin_inset Float figure
wide false
sideways false
status open

\begin_layout Plain Layout
\align center
\begin_inset Graphics
	filename Images/implot-betelgueuse-perpendicular.png
	lyxscale 50
	scale 35

\end_inset


\end_layout

\begin_layout Plain Layout
\begin_inset Caption

\begin_layout Plain Layout
Línia a la meitat de l'espectre (blanc), a l'inici de l'espectre (vermell)
 i al final (verd).
\end_layout

\end_inset


\end_layout

\end_inset

Situar el cursor sobre cada pic i polsant 'espai' s'imprimeix per pantalla
 les coordenades, comparant les X s'obté que el desplaçament del pic es
 de 8 píxels entre l'inici i el final de l'espectre, prou raonable com per
 no haver de fer cap modificació a la imatge.
\end_layout

\end_deeper
\begin_layout Subsection
Extracció de l'espectre
\end_layout

\begin_layout Subsubsection
\begin_inset CommandInset label
LatexCommand label
name "sub:Mode-interactiu"

\end_inset

Mode interactiu
\end_layout

\begin_layout Enumerate
Es carreguen els paquets adients:
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

ecl> imred
\end_layout

\begin_layout Plain Layout

imred> ctioslit
\end_layout

\end_inset


\end_layout

\begin_layout Enumerate
La comanda 'apall' serà l'utilitzada per extreure l'espectre amb els següents
 paràmetres:
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

display grup3_baixa_obj_betelgeuse_30_1.fits 1 fill+
\end_layout

\begin_layout Plain Layout

\end_layout

\begin_layout Plain Layout

apall grup*_*_obj_*.fits line=INDEF nsum=20 width=20 b_sample='-100:-50,50:100'
 b_naver=-100 b_funct=legendre b_order=1 t_nsum=20 t_step=10 t_funct=spline3
 t_order=4 t_niter=1 background=fit weights=none clean=no format=multispec
 nfind=1 find+ recenter- resize- edit+ trace+ fittrac+ extract+ extras+
 review- interactive+
\end_layout

\end_inset


\end_layout

\begin_deeper
\begin_layout Enumerate
Apertura:
\end_layout

\begin_deeper
\begin_layout Enumerate
'line=INDEF': En comptes d'especificar manualment la línia a la que es troba
 l'espectre, IRAF la cercarà automàticament.
\end_layout

\begin_layout Enumerate
'nsum=20': Nombre de línies de dispersió que es sumaran per trobar el centre
 del paràmetre anterior.
 Es recomana que el seu valor coincideixi aproximadament amb la FHWM de
 la traça (extret mitjançant 'implot' tal i com s'especifica a la secció
 
\begin_inset CommandInset ref
LatexCommand ref
reference "sub:Característiques-de-les"

\end_inset

).
\end_layout

\begin_layout Enumerate
'width=20': Tamany del perfil de l'estrella de base a base (igual que l'anterior
).
\end_layout

\begin_layout Enumerate
'resize=yes': Permet que IRAF identifiqui automàticament els límits de la
 traça.
\end_layout

\end_deeper
\begin_layout Enumerate
Finestra de fons de cel:
\end_layout

\begin_deeper
\begin_layout Enumerate
b_sample='-100:-50,50:100': Ubicació del fons de cel relatiu al centre del
 perfil de l'estrella.
 Els valors adients es poden extreure a partir de l'anàlisi presentat a
 la secció 
\begin_inset CommandInset ref
LatexCommand ref
reference "sub:Característiques-de-les"

\end_inset

.
\end_layout

\begin_layout Enumerate
'b_funct=legendre', 'b_order=1' i 'b_naver=-100': Funció utilizada per ajustar
 les mesures del fons de cel.
\end_layout

\end_deeper
\begin_layout Enumerate
Control de la traça:
\end_layout

\begin_deeper
\begin_layout Enumerate
't_nsum=20' i 't_step=10': Nombre de línies que seran utilitzades per identifica
r les variacions de la traça (en cas de que aquesta no sigui perfectament
 recta).
 Es recomana un valor similar al paràmetre 'nsum'.
\end_layout

\begin_layout Enumerate
't_funct=spline3, 't_order=4' i 't_niter=1': Funció que s'utilitzarà per
 ajustar la traça.
\end_layout

\end_deeper
\begin_layout Enumerate
Extracció del fons de cel:
\end_layout

\begin_deeper
\begin_layout Enumerate
'background=fit' i 'weights=none': Restar el fons de cel a la traça.
\end_layout

\end_deeper
\begin_layout Enumerate
Format del fitxer de sortida:
\end_layout

\begin_deeper
\begin_layout Enumerate
'format=multispec' i 'extras=yes': En un mateix fitxer es guarden diverses
 imatges.
\end_layout

\begin_deeper
\begin_layout Enumerate
Banda 1: Espectre de l'estrella.
\end_layout

\begin_layout Enumerate
Banda 2: Espectre del cel.
\end_layout

\end_deeper
\end_deeper
\end_deeper
\begin_layout Enumerate
Al executar la comanda 'apall' del punt anterior, el sistema realitza dos
 preguntes a les que s'ha d'indicar 'yes': 
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

Find apertures for grup3_baixa_obj_betelgeuse_30_1?  (yes): yes
\end_layout

\begin_layout Plain Layout

Edit apertures for grup3_baixa_obj_betelgeuse_30_1?  (yes): yes
\end_layout

\end_inset

A continuació s'obrira l'editor d'apertura, des d'on també es podrà accedir
 a l'editor de l'ajustament de cel.
 Les següents comandes poden ser utilitzades sobre la finestra:
\end_layout

\begin_deeper
\begin_layout Enumerate
A l'editor d'apertura:
\begin_inset Float figure
wide false
sideways false
status open

\begin_layout Plain Layout
\align center
\begin_inset Graphics
	filename Images/Editor apertura.png
	lyxscale 50
	scale 35

\end_inset


\end_layout

\begin_layout Plain Layout
\begin_inset Caption

\begin_layout Plain Layout
Editor d'apertura
\end_layout

\end_inset


\end_layout

\end_inset


\end_layout

\begin_deeper
\begin_layout Enumerate
'y': Canvia l'amplada de l'apertura en funció de la posició del cursor en
 l'eix de les 'y'.
\end_layout

\begin_layout Enumerate
'b': Obre l'editor de l'ajustament de fons de cel.
\end_layout

\begin_layout Enumerate
'q': Finalitza.
\end_layout

\end_deeper
\begin_layout Enumerate
Al l'editor de l'ajustament de fons de cel:
\begin_inset Float figure
wide false
sideways false
status open

\begin_layout Plain Layout
\align center
\begin_inset Graphics
	filename Images/Editor fons de cel.png
	lyxscale 50
	scale 35

\end_inset


\end_layout

\begin_layout Plain Layout
\begin_inset Caption

\begin_layout Plain Layout
Editor de l'ajustament de fons de cel
\end_layout

\end_inset


\end_layout

\end_inset


\end_layout

\begin_deeper
\begin_layout Enumerate
't': Oblida les regions de fons de cel actuals.
\end_layout

\begin_layout Enumerate
's': Dos vegades amb el punter a posicions diferents per tal de marcar noves
 regions a ser considerades fons de cel.
\end_layout

\begin_layout Enumerate
'f': Ajusta el fons de cel amb les regions actuals.
\end_layout

\begin_layout Enumerate
'q': Torna a l'editor d'apertura.
\end_layout

\end_deeper
\end_deeper
\begin_layout Enumerate
Al finalitzar el punt anterior pressionant 'q' a l'editor d'apertura, el
 sistema realitza tres preguntes a les que s'ha de respondre 'yes':
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

Trace apertures for grup3_baixa_obj_betelgeuse_30_1?  (yes): yes
\end_layout

\begin_layout Plain Layout

Fit traced positions for grup3_baixa_obj_betelgeuse_30_1?  (yes): yes
\end_layout

\begin_layout Plain Layout

Fit curve to aperture 1 of grup3_baixa_obj_betelgeuse_30_1?  (yes): yes
\end_layout

\end_inset

A continuació visualitza el centre de la traça i la funció ajustada.
 Aquesta pot ser modificada amb les comandes:
\begin_inset Float figure
wide false
sideways false
status open

\begin_layout Plain Layout
\align center
\begin_inset Graphics
	filename Images/Editor traca.png
	lyxscale 50
	scale 35

\end_inset


\end_layout

\begin_layout Plain Layout
\begin_inset Caption

\begin_layout Plain Layout
Editor de l'ajustament de la traça
\end_layout

\end_inset


\end_layout

\end_inset


\end_layout

\begin_deeper
\begin_layout Enumerate
':spine3', 'legendre': Canviar la funció.
\end_layout

\begin_layout Enumerate
':ordre 2': Canviar l'ordre de la funció.
\end_layout

\begin_layout Enumerate
'f': Ajustar la funció.
\end_layout

\begin_layout Enumerate
'q': Finalitza.
\end_layout

\end_deeper
\begin_layout Enumerate
Finalment, el procés realitza quatre preguntes que s'han de respondre afirmativa
ment: 
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

Write apertures for grup3_baixa_obj_betelgeuse_30_1 to database?  (yes):
 yes
\end_layout

\begin_layout Plain Layout

Extract aperture spectra for grup3_baixa_obj_betelgeuse_30_1?  (yes): yes
\end_layout

\begin_layout Plain Layout

Review extracted spectra from grup3_baixa_obj_betelgeuse_30_1?  (yes): yes
\end_layout

\begin_layout Plain Layout

Review extracted spectrum for aperture 1 from grup3_baixa_obj_betelgeuse_30_1?
  (yes): yes
\end_layout

\end_inset

L'espectre obtingut és visualitzar i per sortir només serà necessari pressionar
 'q':
\begin_inset Float figure
wide false
sideways false
status open

\begin_layout Plain Layout
\align center
\begin_inset Graphics
	filename Images/Visualitza espectre.png
	lyxscale 50
	scale 35

\end_inset


\end_layout

\begin_layout Plain Layout
\begin_inset Caption

\begin_layout Plain Layout
Visualitzador de l'espectre
\end_layout

\end_inset


\end_layout

\end_inset


\end_layout

\begin_layout Enumerate
La informació del procés d'extracció es guarda al directori 'database' amb
 el prefix 'ap' i el nom del fitxer original.
 L'espectre s'ha emmagatzemat al fitxer 'grup3_baixa_obj_betelgeuse_30_1.ms.fits'.
 En cas de voler tornar a visualitzar-ho s'ha d'executar:
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

splot grup3_baixa_obj_betelgeuse_30_1.ms.fits band=1
\end_layout

\end_inset

Si per altra banda, es vol visualitzar l'espectre del fons de cel extret:
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

splot grup3_baixa_obj_betelgeuse_30_1.ms.fits band=2
\end_layout

\end_inset


\end_layout

\begin_layout Subsubsection
Mode automàtic
\end_layout

\begin_layout Standard
En cas de que prèviament s'hagi realitzat un ajustament i extracció interactiva
 de l'espectre, és possible tornar a extreure'l utilitzant els valors guardats
 al directori 'database': 
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

apall grup*_*_obj_*.fits format=multispec recenter- find- resize- edit- trace+
 fittrac- background=fit extract+ extras+ review- interactive-
\end_layout

\end_inset


\end_layout

\begin_layout Subsection
Calibració
\end_layout

\begin_layout Standard
Per tal d'associar una longitud d'ona a cada píxel s'han de realitzar les
 següents accions:
\end_layout

\begin_layout Enumerate
En primer lloc, s'ha d'extreure l'espectre de la lampara amb les mateixes
 condicions que la d'un objecte que s'hagi tractat prèviament segons la
 secció 
\begin_inset CommandInset ref
LatexCommand ref
reference "sub:Mode-interactiu"

\end_inset

:
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

apall grup3_baixa_lamp_10_1.fits out=grup3_baixa_obj_betelgeuse_30_1.ms.lamp.fits
 ref=grup3_baixa_obj_betelgeuse_30_1.fits recenter- trace- background- interactiv
e- nfind=1
\end_layout

\end_inset

Aquest procés crea el fitxer 'grup3_baixa_obj_betelgeuse_30_1.ms.lamp.fits'.
\end_layout

\begin_layout Enumerate
S'utilitzarà un llistat de línies que produeix la lampara, IRAF incorpora
 diversos llistat que poden ser consultats amb:
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

page linelists$README
\end_layout

\end_inset

No obstant, pel cas actual es creara un fitxer al directori 'imdir' anomenat
 'lampara.dat' amb un llistat de longituds (angstroms) i descripcions:
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

# units Angstroms
\end_layout

\begin_layout Plain Layout

4054.0    mercury
\end_layout

\begin_layout Plain Layout

4366.0    mercury
\end_layout

\begin_layout Plain Layout

4877.0    terbium from Tb3+
\end_layout

\begin_layout Plain Layout

5424.0    terbium from Tb3+
\end_layout

\begin_layout Plain Layout

5465.0    mercury
\end_layout

\begin_layout Plain Layout

5777.0    likely terbium from Tb3+ or mercury
\end_layout

\begin_layout Plain Layout

5802.0    mercury or terbium from Tb3+
\end_layout

\begin_layout Plain Layout

5840.0    possibly terbium from Tb3+ or europium
\end_layout

\begin_layout Plain Layout

5876.0    likely europium in Eu+3:Y2O3
\end_layout

\begin_layout Plain Layout

5934.0    likely europium in Eu+3:Y2O3
\end_layout

\begin_layout Plain Layout

5997.0    likely europium in Eu+3:Y2O3
\end_layout

\begin_layout Plain Layout

6116.0    europium in Eu+3:Y2O3
\end_layout

\begin_layout Plain Layout

6257.0    likely terbium from
\end_layout

\begin_layout Plain Layout

6311.0    likely europium in Eu+3:Y2O3
\end_layout

\begin_layout Plain Layout

6508.0    likely europium in Eu+3:Y2O3
\end_layout

\end_inset


\end_layout

\begin_layout Enumerate
A l'espectre estret al primer punt, es necessari identificar associar-ho
 amb com a mínim un parell de línies amb la seva corresponent longitud d'ona:
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

identify grup3_baixa_obj_betelgeuse_30_1.ms.lamp.fits coordlist=imdir$lampara.dat
 fwidth=64
\end_layout

\end_inset


\begin_inset Float figure
wide false
sideways false
status open

\begin_layout Plain Layout
\align center
\begin_inset Graphics
	filename Images/Lamp Calibration.png
	lyxscale 50
	scale 35

\end_inset


\end_layout

\begin_layout Plain Layout
\begin_inset Caption

\begin_layout Plain Layout
Calibració de la lampara.
\end_layout

\end_inset


\end_layout

\end_inset


\end_layout

\begin_deeper
\begin_layout Enumerate
'm': Marcar una línia i introduir a continuació la longitud d'ona a la que
 correspon en Angstroms.
 Per exemple, les més intenses de la lampara utilitzada:
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

977.53 977.52698 (     INDEF ): 6116 
\end_layout

\begin_layout Plain Layout

821.03 821.02771 (     INDEF ): 5465 
\end_layout

\end_inset


\end_layout

\begin_layout Enumerate
'l': Intentar assignar automàticament la resta de línies segons el llistat
 de coordenades especificat.
\end_layout

\begin_layout Enumerate
'f': Mostrar l'ajustament.
\end_layout

\begin_layout Enumerate
'q': Finalitzar i guardar al directori 'database': 
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

Write feature data to the database (yes)? 
\end_layout

\end_inset


\end_layout

\end_deeper
\begin_layout Enumerate
Tots els espectres d'objecte s'han d'associar amb l'espectre de la lampara
 mitjançant un camp de la capçalera:
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

hedit grup3_baixa_obj_betelgeuse_30_1.ms.fits REFSPEC1 "grup3_baixa_obj_betelgeuse
_30_1.ms.lamp.fits" add+ ver- show+
\end_layout

\end_inset


\end_layout

\begin_layout Enumerate
Finalment, es transforma l'espectre de l'objecte per tal de que l'eix d'abscisse
s disposi de la longitud d'ona adient:
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

dispcor grup3_baixa_obj_betelgeuse_30_1.ms.fits grup3_baixa_obj_betelgeuse_30_1.ms.f
its
\end_layout

\end_inset

Per tal de visualitzar el resultat:
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

splot grup3_baixa_obj_betelgeuse_30_1.ms.fits band=1
\end_layout

\end_inset


\end_layout

\begin_layout Enumerate
De forma opcional, és possible normalitzar l'espectre per poder identificar
 més fàcilment les línies d'emissió i/o absorció:
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

continuum grup3_baixa_obj_betelgeuse_30_1.ms.fits output=grup3_baixa_obj_betelgeus
e_30_1.ms.norm.fits function=spline3 order=8 
\end_layout

\end_inset

Algunes de les funcions disponibles a la finestra del gràfic:
\end_layout

\begin_deeper
\begin_layout Enumerate
':function spline3' i ':order 8': Canviar la funció ajustada (cos negre).
\end_layout

\begin_layout Enumerate
'k': Visualitzar espectre normalitzat.
\end_layout

\begin_layout Enumerate
'h': Visualitzar espectre sense normalitzar.
\end_layout

\begin_layout Standard
Per tal de visualitzar el resultat:
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

splot grup3_baixa_obj_betelgeuse_30_1.ms.norm.fits band=1
\end_layout

\end_inset


\end_layout

\end_deeper
\begin_layout Subsection
Mesurar línies
\end_layout

\begin_layout Standard
Un cop obtingut l'espectre, pot ser d'interès analitzar les línies d'absorció
 i d'emissió que apareixen per tal d'identificar la longitud d'ona a la
 que s'han registrat.
\end_layout

\begin_layout Enumerate
Per analitzar un espectre es farà ús de la comanda 'splot':
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

splot grup3_baixa_obj_betelgeuse_30_1.ms.fits band=1 nerrsample=100 sigma0=2.94
 invgain=0.40 save_file="splot.log"
\end_layout

\end_inset

Els paràmetres anteriors indiquen:
\end_layout

\begin_deeper
\begin_layout Enumerate
'nerrsample=100': Activa el càlcul d'errors.
\end_layout

\begin_layout Enumerate
'sigma0=2.94' i 'invgain=0.40': La desviació dels píxels és modelada amb
\begin_inset Formula $\sigma^{2}=\sigma_{0}^{2}+invgain\times I$
\end_inset

 a on 
\begin_inset Formula $I$
\end_inset

 és el valor del píxel.
 Pel cas actual es recomana que 'sigma0' i 'invgain' siguin calculats a
 partir de:
\end_layout

\begin_deeper
\begin_layout Enumerate
\begin_inset Formula $\sigma_{0}^{2}=\frac{readnoise}{\sqrt{width_{aperture}}}=\frac{15}{\sqrt{26}}=2.94$
\end_inset

 a on l'amplada de l'apertura es pot trobar al fitxer 'database/apgrup3_baixa_ob
j_betelgeuse_30_1' ( 
\begin_inset Formula $width_{aperture}\simeq16-\left(-10\right)$
\end_inset

):
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

        low     -764.
 -10.88208
\end_layout

\begin_layout Plain Layout

        high    765.
 16.75909
\end_layout

\end_inset


\end_layout

\begin_layout Enumerate
\begin_inset Formula $invgain=\frac{1}{gain}=\frac{1}{2.5}=0.40$
\end_inset

 a on el guany ve determinat per les característiques de la CCD (secció
 
\begin_inset CommandInset ref
LatexCommand ref
reference "sec:Instrumentació"

\end_inset

).
\end_layout

\end_deeper
\begin_layout Standard
\begin_inset Float figure
wide false
sideways false
status open

\begin_layout Plain Layout
\align center
\begin_inset Graphics
	filename Images/Analisi splot.png
	lyxscale 50
	scale 35

\end_inset


\end_layout

\begin_layout Plain Layout
\begin_inset Caption

\begin_layout Plain Layout
Anàlisi d'espectres amb 'splot'.
\end_layout

\end_inset


\end_layout

\end_inset

A la finestra de la gràfica es poden utilitzar les següents funcions:
\end_layout

\begin_layout Enumerate
'a': Polsant 2 cops realitza un zoom entre les dos ubicacions a on estigués
 el punter.
\end_layout

\begin_layout Enumerate
'c': Torna a visualitzar l'espectre anul·lant el zoom.
\end_layout

\begin_layout Enumerate
'espai': Imprimeix per pantalla la posició actual del punter.
\end_layout

\begin_layout Enumerate
'k': Polsant 2 vegades ajusta una gausiana simple.
 El punter ha d'estar per sota si és una línia d'emissió i per sobre si
 és d'absorció.
 En cas que es vulguin ajustar regions amb diferents pic és recomana l'ús
 de 'd'.
\end_layout

\begin_layout Enumerate
'd': Polsant 2 vegades es selecciona la regió que es vol ajustar a una de
 les següents funcions: 
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

Lines ('f'ile, 'g'aussian, 'l'orentzian, 'v'oigt, 't'ype, 'q'uit:
\end_layout

\end_inset

Per exemple, en el cas de les gaussianes, s'haurà de pressionar 'g' situant
 el cursor sobre cadascun dels pics de la gràfica.
 Un cop finalitzat, es pressiona 'q' i es respon 'a', 'a' i 'y' a les preguntes:
 
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

Fit positions (fixed, single, all, quit)
\end_layout

\begin_layout Plain Layout

Fit Gaussian widths (fixed, single, all, quit)
\end_layout

\begin_layout Plain Layout

Fit background (no, yes, quit)
\end_layout

\end_inset

Sota del gràfic es mostra informació sobre cada pic senyalat, es pot rotar
 entre ells amb les tecles '+' i '-'.
 Cal pressionar 'q' dos vegades per sortir i 'c' per esborrar les línies
 ajustades.
\end_layout

\begin_layout Enumerate
'e': Polsant 2 vegades es selecciona la regió sobre la qual es calcula el
 punt mig, entre altres paràmetres.
\end_layout

\begin_layout Enumerate
'm': Polsant 2 vegades calcula la mitja, RMS i relació senyal soroll (S/N)
 d'una regió.
\end_layout

\begin_layout Enumerate
'$': Intercanvia coordenades físiques (longitud d'ona) i de la imatge (píxels)
\end_layout

\begin_layout Enumerate
':log' i ':nolog': Activa/desactiva el log.
\end_layout

\begin_layout Enumerate
'# text': Afegeix un comentari al log.
\end_layout

\begin_layout Enumerate
':label text', ':mabove text' i ':mbelow text': Afegeix text a on esta situat
 el punter en sentit horitzontal o vertical amb una marca.
\end_layout

\begin_layout Enumerate
'q': Finalitza l'anàlisi.
\end_layout

\end_deeper
\begin_layout Enumerate
Els resultats dels ajustaments realitzat amb 'splot' queden registrats al
 fitxer 
\begin_inset Quotes eld
\end_inset

splot.log
\begin_inset Quotes erd
\end_inset

:
\begin_inset listings
lstparams "basicstyle={\scriptsize},breaklines=true"
inline false
status open

\begin_layout Plain Layout

Jun  5  2:03 [grup3_baixa_obj_betelgeuse_30_1.ms.fits[*,1,1]]:         
\end_layout

\begin_layout Plain Layout

    center      cont      flux       eqw      core     gfwhm     lfwhm
\end_layout

\begin_layout Plain Layout

  5868.464  11581.12  -793614.
     68.53  -4566.84     163.3        0.
\end_layout

\begin_layout Plain Layout

  (0.12939) ( 5.9505) (  2948.) ( 0.2343) (  10.74) ( 0.7415) (     0.)
\end_layout

\begin_layout Plain Layout

  6179.632  14239.82  -774879.
     54.42  -5270.39     138.1        0.
\end_layout

\begin_layout Plain Layout

  (0.11816) ( 6.1139) (  2890.) ( 0.1987) (  9.852) ( 0.4916) (     0.)
\end_layout

\begin_layout Plain Layout

    center      cont      flux       eqw      core     gfwhm     lfwhm
\end_layout

\begin_layout Plain Layout

  7418.592  5027.771  3312291.
    -658.8   9111.94     341.5        0.
\end_layout

\begin_layout Plain Layout

  (0.17969) (       ) (  4824.) ( 0.9594) (  6.614) ( 0.5933) (     0.)
\end_layout

\end_inset


\end_layout

\end_body
\end_document
